Origin of the photoemission final-state effects in Bi 2 Sr 2 CaCu 2 08 
by very-low-energy electron diffraction 
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Very-low-energy electron diffraction with a support of full- 
potential band calculations is used to achieve the energy po- 
sitions, K|| dispersions, lifetimes and Fourier compositions of 
the photoemission final states in Bi2Sr2CaCu2 0s at low ex- 
citation energies. Highly structured final states explain the 
dramatic matrix element effects in photoemission. Intense 
c(2x2) diffraction reveals a significant extrinsic contribution 
to the shadow Fermi surface. The final-state diffraction ef- 
fects can be utilized to tune the photoemission experiment on 
specific valence states or Fermi surface replicas. 



Photoemission (PE) spectra of Bi 2 Sr 2 CaCu208 
(Bi2212) and other high-temperature superconductors 
can significantly deviate from the photohole spectral 
function A(k, w) due to an energy and k dependence of 
the matrix element (ME) [1-4]. For example, the long 
debated dependence of the shape of the observed Fermi 
surface (FS) of Bi2212 on the photon energy hv and ex- 
perimental geometry has finally found its explanation in 
ME effects. A predictive use of these effects can help to 
resolve subtle details of the electronic structure of Bi2212 
such as the bilayer splitting [2,3]. Presently, analysis of 
ME effects employs computational simulations of the PE 
process [1,2]. However, there remains an inherent uncer- 
tainty of a few eV due to the lack of any reliable knowl- 
edge about the energy location and lifetime of the PE 
final states. Additional complications are incommensu- 
rate (IC) structural modulations and antiferromagnetic 
(AFM) fluctuations, which manifest themselves in PE ex- 
periments, respectively, as diffraction replicas of the FS 
and so-called shadow FSs [5-7] . 

Very-Low-Energy Electron Diffraction (VLEED) gives 
direct access to the PE final states (sec [8,9] and the ref- 
erences therein) based on the fact that, neglecting the 
electron-hole interaction within the sudden approxima- 
tion, the PE final states are the time-reversed LEED 
states [10]. In the VLEED elastic electron reflectivity 
spectra R{E) , the energies of the spectral structures give 
the final-state dispersions E(k), and their broadening 
the corresponding lifetimes r expressed by the imaginary 
part ImS of the self-energy as H/t = 2ImE. 

We here use VLEED to achieve first direct information 



on the PE final states in Bi2212 up to 45 eV above the 
Fermi level Ep and elucidate its immediate implications 
for the PE experiment. 

Experiment and band calculations - The VLEED ex- 
periment was carried out using a standard LEED setup 
operated in the retarding field mode [9,11]. The total 
R(E) integrated over all diffracted beams (plus structure- 
less inelastic reflectivity) was measured in the target cur- 
rent T(E). Intensities of individual beams were measured 
on the LEED screen using a CCD camera. Three sam- 
ples with an oxygen content near optimal doping showed 
identical results. 

The band calculations employed the full-potential 
APW+local orbitals method as implemented in the 
WIEN2k package [12]. The generalized gradient approxi- 
mation (GGA) was used to describe exchange-correlation 
within the density functional theory (DFT). The calcula- 
tions employed the body-centered-tetragonal (bet) crys- 
tal structure and a basis of 2200 augmented plane waves. 
The obtained E(k), characterizing essentially the bulk 
states with infinite lifetime, was connected to the VLEED 
process based on Fourier expansion of the Bloch waves, 
which was used to evaluate the partial currents Tk into 
the crystal excited in each Bloch wave, and the integral 
T(E) spectra [15]. Finite electron lifetime was simulated 
by Lorentzian smoothing and damping of T(E) [8]. 

The experimental structural parameters from [13] and 
[14] were tested in the calculations, with the latters yield- 
ing the best agreement with the VLEED experiment. 
While the valence bands changed insignificantly, the un- 
occupied E(k) was strongly sensitive to the structural 
parameters. Such a sensitivity suggests that any model 
of the final states, used in PE simulations, needs to be 
checked against the VLEED experiment. 

Integral T(E) spectra: Dispersions and lifetimes of the 
final states - The structures in T(E) all reflect those in 
elastic R(E), with the inelastic reflectivity giving only 
a featureless background. We start with analysis of the 
normal incidence spectrum in Fig. 1(a). It corresponds, 
by parallel momentum conservation and assuming in- 
significant umklapp contributions from the IC super- 
structure, to the TZ direction of the Brillouin zone (BZ). 
A slope of T(E) towards higher energies is due to an in- 
crease of the inelastic reflection. Its unusual steepness 
may indicate strong electron coupling in the cascade re- 
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laxation process. To emphasize the T(E) structures, we 
subtracted a quadratic background, Fig.l(&). 

The calculated T(E) spectrum (minus a quadratic 
background) is also shown in Fig. 1(6). A reasonable 
agreement with the experiment, especially considering 
the neglect of umklapp, confirms our theoretical model 
of the unoccupied states. The best averall agreement is 
achieved by shifting the calculations by ~3 eV to higher 
energies. Besides the computational approximations, this 
can be attributed to excited-state self-energy corrections 
AE to the DFT-GGA exchange-correlation. The experi- 
ment suggests certain band dependence of AE. 

Guided by the calculated T(E), we can trace the 
band structure origin of the experimental T(E) struc- 
tures. The T(E) maxima are seen to reflect manifolds 
of so-called coupling bands (= dominant final states in 
PE) which arc characterized by large partial Tk currents 
[9,15]. The T(E) minima reflect the gaps between such 
manifolds. Individual bands arc however not resolved in 
T(E) because their energy separation is much smaller 
than ImE. Note that the coupling bands, apart from ef- 
fective coupling to the incident plane wave, must have a 
3-dimcnsional (3D) character with considerable k± dis- 
persion to transmit current into the crystal; they origi- 
nate from the electronic orbitals oriented off-plane. The 
2-dimcnsional (2D) bands, originating from the in-plane 
orbitals, arc ineffective in both VLEED and PE. 

Finite electron lifetime manifests itself in broaden- 
ing and damping of the T(E) structures. We esti- 
mated the energy dependence of Im E based on its ap- 
proximate parametric relation to the dielectric function 
e as ImS(£) = a + b Im [— 1/e (q = 0, u)] dui with 
e (q = 0, u) taken from the electron energy loss data [16]. 
The parameters a and b were optimized to bring the cal- 
culated T{E) to the best agreement with the experiment 
in broadening and damping of the spectral structures. 
The obtained ImE(£ l ), which yields the calculated T(E) 
in Fig.l (6), is shown in the insert in Fig.l (a) super- 
imposed on e(q = 0,w). The increase of ImE results 
from two plasmon peaks in e (q = 0,w) associated with 
the CuO valence bands [16]. This increase of ImE in the 
final states explains, for example, why the PE intensity 
ratio of the CuO bilayer bonding and antibonding states 
varies with hv above 30 eV much smoother in the exper- 
iment than in simulations employing constant ImE =1 
eV [17] relevant closer to hv of 20 eV. 

Based on the experimental ImE(£) and calculated 
band dispersions, we found that the inelastic electron 
scattering limits the mean free path A around 3A, de- 
creasing with energy, in agreement with the electron en- 
ergy loss data [16]. In the final-state band gaps A is fur- 
ther reduced by elastic scattering off the crystal potential, 
resulting in a diffraction structure in X(E), particularly 
strong where ImE is small. Qualitatively, it follows the 
structure in T(E) where the maxima reveal the 3D states 



enabling effective electron transport (accurate evaluation 
of A involves calculations of damped Bloch waves [9,11]). 
Predictive tuning A based on the experimental T(E) can 
therefore be employed in the PE experiment to resolve 
signals from different atomic layers. 

Experimental dispersion of the T(E) spectra with the 
incidence Km varying in the TX azimuth (perpendicular 
to the IC superstructure) is shown in Fig. 2 (left). The 
shading shows the energy intervals of the T(E) maxima 
reflecting the coupling bands. This map gives therefore 
the surface-projected K|| dispersions of the final states. 

Our VLEED results give the first direct evidence of 
unoccupied quasiparticle dispersions in Bi2212, whereas 
the inverse PE studies near Ep (see, e.g., [18]) found 
only dispersionless states. Furthermore, the highly struc- 
tured non-parabolic dispersions illustrate a general fact 
that at low excitation energies the final states, especially 
for complex materials, dramatically deviate from a free- 
electron-like continuum due to multiple scattering on the 
crystal potential. This property of the final states is the 
essential origin of dramatic hv dependences in PE. 

The results for the TY azimuth (along the IC super- 
structure) are shown in Fig. 2 (right). By the fundamental 
lattice symmetry the TY azimuth is equivalent to TX, 
but this is complicated by the IC superstructure. The 
ry dispersion image can then be viewed as the TX dis- 
persion with its IC umklapp replicas superimposed and 
hybridized with the main dispersion. The observed image 
resembles therefore a smeared TX dispersion. 

The dispersion in TY comes at first glance as a sur- 
prise: Competing fundamental and IC potentials results 
in the collapse of strict periodicity (and the entire band 
structure concept) and thereby the parallel momentum 
conservation. As all states are now available for any in- 
cident K||, one might expect that any dispersion in K|| 
would be lost. However, the VLEED spectra depend, 
apart from bare availability of the states, on how their 
wavefunctions match the incident plane wave [15]. Vari- 
ation of the matching conditions with results in the 
observed dispersion. A similar appearance of competing 
periodicities in PE was elucidated in [19]. 

Our experimental T(E) data yield thus the energy 
positions, K|j dispersions and lifetimes of the PE final 
states. Apart from direct implications in the PE experi- 
ment such as tuning A in energy and K|| , this information 
can be used as a vital input in computational simulations 
of the ME effects instrumental for resolving details of the 
valence E(\s.) [2]. ME changes are expected between the 
T(E) maxima owing to a phase change in the final-state 
wavefunction through the band gaps. 

Intensities of diffracted beams: hidden c(2x2) periodic- 
ity and final-state diffraction effects - A typical normal- 
incidence VLEED image is shown in Fig. 3 (left). Due to 
the IC superstructure the diffracted intensity is arranged 
in lines along TY [18,20], with non- vanishing intensity 
between the spots manifesting the fractal character of 
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diffraction in presence of two competing periodicities. 
The intense spots, each with IC satellites, are seen to 
reflect, firstly, the reciprocal lattice of the fundamental 
lattice at (±27r, 0) and (0, ±2w). Surprisingly, there are 
other intense spots (±7r, ±7r) (arrows), although weaker 
in (— it, it) and Their symmetry reveals a hid- 

den c(2x2) periodicity This phenomenon is observed in 
electron diffraction for the first time. 

I(V) data are shown in Fig.3 {right) as the VLEED in- 
tensity profile along an IC diffraction line, as a function 
of energy and qi-K,-*) momentum transfer. All diffrac- 
tion spots show a significant and non-monotonous energy 
dependence beyond the kinematic picture. The c(2x2) 
spots are intense only below ~20 eV, and in the normal 
LEED energy range they vanish. 

To identify the origin of the c(2x2) periodicity, we first 
performed a kinematic simulation using the IC super- 
structure parameters from [21]. It has however yielded 
only the fundamental-lattice spots with IC satellites. Be- 
yond the kinematic picture, possible origins of the c(2x2) 
spots should be similar to the shadow FS with its c(2x2) 
periodicity: either two inequivalent Cu sites per plane 
in the full orthorhombic unit cell of Bi2212 [6], or short- 
range AFM correlations within the CuO planes [5] mani- 
festing themselves in VLEED through the exchange scat- 
tering [22]. However, the observed K|| broadening of both 
c(2x2) spots and shadow FS seems too small compared 
to the short AFM correlation lengths found by neutron 
scattering. Finally, the c(2x2) spots can be peaks of mul- 
tiple scattering on the competing fundamental and IC 
potentials. Whichever the origin, the diffraction process 
involves multiple scattering between the CuO planes, and 
the outer SrO and BiO planes with their IC modulations 
[21]. The involvement of the CuO planes is corroborated 
the c(2x2) intensity reduction with energy which goes 
along with the decrease of A. 

Complementing the information on the final states ob- 
tained from T(E), the VLEED intensity distribution re- 
flects their Fourier composition: Being the time-reversed 
LEED state, the final state can be expanded over the 
surface reciprocal vectors $ (r) = <fi % (r±) e*( k||+s ) r|1 , 
where the Fourier coefficients <fi % (r±) correspond, in vac- 
uum, to the incident and diffracted plane waves and, in 
the crystal, to the Bloch wavefield. The VLEED inten- 
sities Ig are the vacuum asymptotics of \(j> s (r±)\ 2 . As 
(fi s (r± ) run across the surface continuously in ampli- 
tude and derivative, I s reflect, approximately, the Fourier 
composition of the final state in the crystal. 

In the PE experiment this information immediately 
characterizes the final-state diffraction effects. For exam- 
ple, in the FS mapping experiment a g-diffraction replica 
of the main FS occurs whenever upon variation of k|| the 
corresponing ky + g Fourier component of the final state 
matches, in the parallel momentum, any Fourier compo- 
nent of the initial state, giving rise to non-zero overlap of 



the final- and initial-state wavefunctions. The PE signal 
is then proportional, apart from phases of the Fourier 
components and selection rules entering the ME [1,2,7], 
to the final-state |<^ g (rj_)| 2 and thus to its 7 g asymp- 
totics. Therefore, the PE signal of certain g-diffraction 
replica follows, roughly, the g-diffracted intensity ob- 
served in the VLEED image. The I(V) data can then 
be used to tune the final-state energy to enhance or sup- 
press certain diffraction replicas in the FS mapping ex- 
periment. Similar ideas can be exploited to single certain 
bands out of the valence band multitude based on their 
dominant Fourier component. Apart from Bi2212, such 
an approach can be applied to detail electronic structures 
of other complex materials, for example, perovskites. 

Our VLEED data on the final-state diffraction effects 
leads to significant revisions in the current interpretation 
of some PE data on Bi2212. For example, the shadow 
FS with its c(2x2) periodicity is commonly considered as 
an intrinsic effect pertinent to the initial state. To check 
against the extrinsic effects due to the final-state diffrac- 
tion, we have measured VLEED images with incidence 
K|| set near the M point where a large shadow FS signal 
is found in the PE experiment. The c(2x2) spots were 
found intense below ~25 eV and gradually vanishing at 
higher energies. This indicates that at low hv, commonly 
used in the PE experiment, the shadow FS intensity has 
a large extrinsic contribution. This is corroborated by 
the absence of a hybridization gap between the main and 
shadow FS [7] . The intrinsic shadow FS signal should be 
measured with higher hv above 25 eV. 

Finally, it should be noted that the final-state energies 
in PE, strictly speaking, differ from those observed in 
VLEED by the electron-hole interaction (excitonic) en- 
ergy E e -h- This figure, characterizing deviations from 
the sudden approximation, can be determined by com- 
paring characteristic points in the PE intensity behavior 
with the predictions of VLEED. For typical layered ma- 
terials such an analysis has however yielded negligible 
E e -h values, at least on the scale of 100 meV [8]. 

Conclusion - Using VLEED with a support of full- 
potential band calculations, we have achieved first direct 
information on the energy positions, K|| dispersions, life- 
times and Fourier composition of the PE final states in 
Bi2212 at the excitation energies up to 45 eV above Ep. 
Implications of this information for the PE experiment 
extend from tuning A to predictive use of the ME ef- 
fects to focus on specific valence states or FS replicas. In 
particular, intense final-state c(2x2) diffraction identifies 
significant extrinsic contribution to the shadow FS. 

We thank M. Golden, A. Zakharov and E.E. Krasovskii 
for promoting discussions. The work was supported by 
Deutsche Forschungsgemeinschaft (grant CI 124/5-1). 
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FIG. 1. (a) Experimental T(E) spectrum; (6) Structure 
in the experimental and calculated T(E) enhanced by sub- 
tracting a quadratic background. The calculation is shifted 
by +3 eV to mimic AE. The inset shows the experimental 
ImS (with the uncertainty bars) compared with Im [— 1/e]; 
(c) Calculated E(k) with the partial currents Tk shown by 
grayscale. The T(E) maxima reflect manifolds of the final 
bands coupling to vacuum. 

FIG. 2. Experimental Ky dispersion of the T(E) spectra. 
The shading shows the energy regions of the T(E) maxima, 
identified by d 2 T/dE 2 > 0. The values of d 2 T/dE 2 in these 
regions are represented in a logarithmic grayscale (black = 
maximal value). This map is the surface-projected dispersion 
of the PE final bands. 

FIG. 3. (left) Normal-incidence VLEED image at an energy 
of 16.4 eV above Ef- The specular beam is in the center ob- 
scured by the electron gun. The diffracted intensity maxima 
reflect the fundamental lattice and a hidden c(2x2) periodic- 
ity (arrows); (right) I(V) data for the indicated IC diffraction 
line shown in a linear grayscale (white = maximal intensity). 
The coordinates on the LEED screen, somewhat distorted by 
the retarding field, are rendered into the q^,-^) values. 
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Fig.l by Strocov et al 
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